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ABSTRACT 


The  acoustic  radiation  from  an  elastic  cylindrical  shell  defined  by 
nonconcentric  cylindrical  surfaces  and  submerged  in  an  acoustic  fluid  is 
studied.  The  driving  mechanism  is  a time  harmonic  internal  pressure.  In 
particular,  the  influence  of  the  nonconcentricity  parameter,  6,  on  the 
resulting  radiation  field  is  examined.  The  solution  is  obtained  by  means 
of  an  expansion  in  powers  of  this  parameter  about  the  concentric  solution, 
6=0.  Each  order  of  solution  is  solved  in  closed  form  by  means  of 
boundary  integral  equation  methods.  Results  are  given  for  the  zeroth  and 
first  order  solutions  and  it  is  seen  that  a driving  mode  will  introduce 
nearest  neighbor  ( (j-1)  and  (j+1)  ) modes  in  the  first  order  solution. 
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I.  INTRODUCTION 

We  consider  the  radiation  into  an  infinite  acoustic  fluid  from  a submerged 
elastic  shell  subjected  to  a time  harmonic  internal  pressure  with  some  pre- 
scribed spatial  distribution.  The  shell  is  bounded  by  surfaces  of  separable 
geometry,  i.e.,  such  that  the  governing  (Helmholtz)  partial  differential  equa- 
tions have  modal  solutions,  that  are  nonconcentric , i.e.,  the  origin  of  the 
coordinate  system  in  which  the  outer  boundary  is  "separable"  is  displaced  from 
that  for  the  inner  boundary  by  a nonconcentricity  distance,  6 . If  6 were  zero, 
the  entire  problem  could  be  solved  by  separation  of  variables  but  with  6 
nonzero,  the  problem  is  not  separable. 

The  purpose  of  this  investigation  is  to  study  the  influence  of  the  noncon- 
centricity on  the  resulting  radiation  field  by  making  an  expansion  of  the 
solution  in  powers  of  6 . Clearly  the  zeroth  order  will  give  the  classical 

separation  of  variables  solution.  Higher  order  solutions  will  then  represent 
the  influence  of  nonconcentricity  on  this  basic  solution. 

Rather  than  work  with  the  governing  differential  equations,  we  convert  to 
a boundary  integral  equation  representation  which  eventually  allows  us  to  solve 
a set  of  algebraic  equations  for  each  mode  in  each  order  to  obtain  the  exact 
solution . 

To  simplify  the  algebra,  the  submerged  elastic  body  will  be  taken  as  an 
infinite  cylinder  in  plane  strain  thereby  reducing  the  problem  to  a two  dimen- 
sional case.  A specific  example  of  nonconcentric  circular  boundaries  is  calcu- 
lated as  an  illustration.  Analogous  three  dimensional  problems  are  also  being 
studied  - the  axisymmetric  three  dimensional  case  requires  essentially  no  more 
effort  than  this  two  dimensional  problem. 

The  full  equations  of  elasticity  shall  be  used  rather  than  approximating 
set  of  shell  equations  since  the  use  of  the  boundary  integral  equation  method 
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reduces  the  problem  by  one  dimension  anyway  thereby  accomplishing  one  of  the  main 
purposes  of  shell  theory. 

This  type  of  problem  has  been  studied  previously  for  acoustic  radiators, 
both  by  the  boundary  integral  equation  method  [1]  and  by  shifting  and  super- 
position theorems  [2]. 


A 
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11.  GENERAL  FORMULATION 

A.  ACOUSTIC  FLUID  EQUATIONS 

The  Infinite  acoustic  fluid  exterior  to  the  elastic  body  is  described 
in  terms  of  an  acoustic  (excess)  pressure  field,  p , which  satisfies  a 
Helmholtz  equation 


(1) 


^7  * ^ + W ^ = o 


where  k = w/c  is  the  acoustic  wave  number,  w is  the  frequency  of  the  time 
harmonic  dependence  (taken  in  the  form  exp[-iiot])  and  c is  the  acoustic  sound 
speed.  This  equation  is  readily  replaced  by  a Helmholtz  (or  Weber,  for  the  two 
dimensional  case)  integral  equation;  e.g.,  [3] 


(2) 


where  G = ■[  H (kR)  in  two  dimensions  and  G = -7  exp[ikRj/R  in  three 

4 o 4 

dimensions.  R is  the  distance  between  the  field  point  r and  the  integration 

variable,  r , n is  the  outward  normal  from  the  fluid,  T is  the  line  or 
o o 

surface  bounding  the  fluid,  i.e.,  the  elastic  solid-acoustic  fluid  interface  with 
length  or  area  element  dS  , and  £ = 0,  1/2,  1 depending  on  whether  the  field 
point  is  exterior  to  the  fluid,  on  the  interface  or  interior  to  the  fluid.  The 
radiation  condition  is  automatically  satisfied  by  the  choice  of  G . Since 
only  Hankel  functions  of  the  first  kind  will  occur,  the  superscript  (1)  will  be 
dropped . 

The  fluid  velocity  field  is  related  to  the  pressure  field  through  the 
acoustic  velocity  potential  41  , which  could  also  have  been  used  as  the  depen- 
dent variable 


<■0  AT 


= - V 4> 


(3) 


(b)  Ip  - jV  'db/dt  ~ -KjJfr  $ 


where  p Is  the  acoustic  fluid  density. 
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B.  K1ASTIC  SHKLL  EQUATIONS 

The  elastic  shell  involves  a finite  domain  with  an  inner  and  an  outer 
boundary  (treating  the  infinite  cylinder  as  a two  dimensional  plane  strain 
problem).  It  is  best  described  by  the  displacement  potentials  <t>  and  S' 
defined  by  the  displacement  field  u as; 

w u.  = V + V*  Y , V-Y-O 

These  potentials  also  satisfy  Helmholtz  equations 

(a)  Si'  + Wo  $ = O 

(b)  V T * Ur  T = O 

where  k = ui/c  , k = a>/c  and  c = J (A+2;i)/p  is  the  dilatational  wave  speed, 

D D K K D S 

c„  = / p/p  is  the  rotational  wave  speed  in  an  elastic  medium  of  density  p 

K s s 

and  Lame  parameters,  A and  p , e.g.,  [A]. 

In  the  two  dimensional  plane  strain  case  considered  here,  the  vector  dis- 
placement potential  S'  has  only  one  nonzero  component  and  will  be  written  as  a 
scalar,  ip  . 

Tlie  displacement  potentials  are  also  related  to  the  dilatation  A and  the 
rotation  lo 

R 

(a)  /\  - - Wo 

(6) 

(b)  CjOr  = * X W 'Y 

which  may  also  be  used  as  dependent  variables. 

Clearly  these  displacement  potentials  also  satisfy  Helmholtz  integral  equa- 
tions analogous  to  equation  (2)  with  two  major  changes.  The  wavenumber  k must 

be  changed  to  either  k or  k and  the  surface  V now  consists  of  two  parts, 

D R 

the  inner  and  the  outer  boundaries.  We  note  that  the  direction  of  n changes 
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sign  between  these  boundaries  since  it  is  required  to  be  "outward"  from  the 
elastic  medium.  Then 


(7)  £§(?)-- 


Po 


d&o 

ar\o 


3 §Cron  \ o' 

2>Ylo  J 0 /JoUTER 


(8) 


* 5 AS...™ 

£ Ycf  )=  § [-  YdO  V??.  + G«  i^r]  A5oUtER 
r0 


t\ 

where  Gp,  GR  are  equivalent  to  the  previous  G with  kp,  kR  replacing 
r in  any  of  these  boundary  integral  equations  will  be  placed  on  (or 
infinitesimally  near  to)  either  the  inner  or  the  outer  surface.  Thus  each 
the  above  boundary  integral  equations  will  produce  two  separate  integral 
equations . 


k . 


of 
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C.  BOUNDARY  CONDITIONS 

We  supplement  these  equations  with  five  boundary  conditions  requiring  zero 
shear  stress,  continuity  of  the  normal  displacement  or  velocity  field  and 
continuity  of  the  normal  stress  field  at  the  outer  boundary  of  the  elastic  shell 
and  zero  shear  stress  and  continuity  of  normal  stress  at  the  inner  boundary. 
(Actually,  those  at  the  outer  boundary  might  better  be  called  interface  condi- 
tions . ) 

For  plane  strain  in  a two  dimensional  orthogonal  curvilinear  coordinate 
system  (n,B),  we  have  the  infinitesimal  strains,  e.g.,  [5] 

r ( 7ft.,)  d + ( ' /PuftO  "U-p  ^ j3 

(9)  (b)  £ pp  = ( '/hi)  + OAuPu)  lie*  3 

<c)  6 o < p = ( fti/d.}  ^ [ Up/dr]/^  + (ft'AO  ^ [U«v/A,]/9p 


and  the  dilatation  and  rotation 


(a)  A = 

(10) 

<b>  GO*  = 

where  h^  and  h 
coordinates,  (fi^ 


[ ^ ( \1<*  0(  + c>  ( Up,  &-i  )/ ^ /[Pud-u] 

[ c>  ( Up9»i)/dof  - 3 (Utf  9l')/  dx] 

? are  the  metric  coefficients  (scale  factors)  for  these 

= 1). 


We  can  also  write  out  equation  (4)  specifically: 

(a)  u*  -- (7M  * C/MaT/dp 

(ID 

(b)  lip  = ('/dx) 

We  shall  choose  the  boundaries  to  be  surfaces  of  constant  a ; the  a 
direction  is  then  normal  to  these  surfaces.  Continuity  of  normal  stress 


requires 
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- = OZ(o<  - /\  A + 2 M- 

which  reduced  to 

„ , t , c afc*  a & * X a*ii  il 

(12)  -ip=-fsw  $ + X5  15  fc?  V sp 

1 ^ f,  i 31.  J V . ± si  3J  4 i t2!f  1 

i hi  a p1  " sv.  M dm  ap  ap  ;><*  w,0u  saapi 


Similarly,  vanishing  of  the  shear  stress  requires 

z ^ $ z 3>fcj  Bj  JL  BPu 

<13)  0 ■'  i ^ ^ ^ a*  3 fa 


_ _L_  I 

^ » 5,  x a oi  a ^ s ^ r* 


and  continuity  of  normal  velocity  requires 


(14)  AJl 
d<* 


f,00l 


1 ^ + j_  3V1 

ft.,  d*  J?Lt  SP  J 
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n SOLUTION  I’ROCFPURK 

There  will  in  general  be  ten  equations,  two  boundary  conditions  plus  two 
integral  equations  (on  T , ij;)  evaluated  at  the  inner  boundary  and  three  boundary 
conditions  plus  three  integral  equations  (on  p,  <!>  and  ilO  evaluated  at  the 
outer  boundary.  Treating  p at  the  inner  boundary  as  known,  these  equations 
describe  the  behavior  of  3i>/3a,  <i»  , 3i|»/3a  on  the  inner  boundary  and 

p,  !p/  iu,  T , 3 4>  / 3 u , ip,  3\p/da  on  the  outer  boundary  as  ten  unknown  functions  of 
the  remaining  independent  variable  B . To  emphasize  the  use  of  a different 
origin,  we  shall  use  the  symbols  (a,B)  for  the  inner  boundary  coordinate 
system  and  the  symbols  (ot',6')  for  the  outer.  We  also  use  a subscript  o to 
represent  an  integration  coordinate.  To  simplify  the  algebra  as  much  as  possible, 
the  inner  boundary  forcing  pressure  will  be  taken  proportional  to  the  jth  eigen- 
function of  the  separated  Helmholtz  equation  for  that  surface.  Finally  the  inner 
boundary  will  be  a surface  of  constant  a = a and  the  outer  of  constant  a'  = 
b,  as  shown  in  Figure  1. 

We  recognize  that  the  inner  boundary  variables  may  be  expanded  as  a series 
in  terms  of  the  inner  boundary  eigenfunctions;  similarly  for  the  outer  boundary 
variables.  The  il  and  3t|)/3«  series  will  be  in  terms  of  the  antisymmetric  modes 
if  p,  l1,  etc.  are  in  terms  of  the  symmetric  modes  and  vice-versa  as  required 
by  the  boundary  conditions.  Since  each  boundary  condition  involves  only  a 
single  boundary,  they  must  be  separable,  i.e. , each  mode  of  the  series 
expansion  must  satisfy  the  boundary  conditions  independently  of  the  other  modes 
thereby  providing  five  algebraic  equations  on  the  unknown  coefficients  of  each 
mode  of  the  expansion. 

The  boundary  integral  equations,  however,  have  one  set  of  terms  wherein  the 
integration  and  field  point  lie  on  the  same  surface,  thereby  allowing  a direct 
integration  to  be  carried  out  independently  of  iS  , and  another  set  of  "mixed" 
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terms  where  the  Integration  and  field  point  lie  on  different  surfaces,  thereby 
involving  5 explicitly. 

We  can  carry  out  the  integration  of  the  first  type  quite  simply  since  the 
Greens'  function  possesses  an  eigenfunction  expansion  in  terms  of  the  same 
functions  (either  completely  inner  or  completely  outer  boundaries)  as  did  the 
other  dependent  variables.  These  integrals  then  pose  no  real  difficulty. 
Integrals  of  the  mixed  (second)  type  would  pose  no  difficulty  if  <$  were  zero, 
i.e.,  in  the  concentric  cases,  the  eigenfunctions  in  terms  of  6 must  be 
identical  on  surfaces  of  constant  a . This  then  suggests  a power  series  expan- 
sion in  terms  of  6 about  6 equal  to  zero  of  all  of  the  algebraic  coefficients 
of  all  of  the  dependent  variable  series  expansions  as  well  as  that  of  the  Greens' 
function  for  the  mixed  integrals.  The  equations  for  the  boundary  conditions  and 
for  the  first  type  of  integral  do  not  involve  6 explicitly  and  expressions  for 
each  order  of  5 are  identical  to  the  original  expression.  The  mixed  integrals 
involve  a product  of  two  series  in  powers  of  6 and  thereby  introduce  a coupl- 
ing of  orders.  The  zeroth  order  solution  will,  of  course,  be  that  obtainable 
directly  by  separation  of  variables  for  the  concentric  case.  The  first  order 
solution  is  coupled  to  the  zeroth  order  solution  through  the  mixed  integrals  — 
this  in  the  first  order  solution  will  be  found  to  introduce  "nearest  neighbor 
modes"  to  each  mode  of  the  zeroth  order  solution,  etc. 

By  taking  the  forcing  function  to  excite  only  a single  mode  in  the  zeroth 
order  solution,  we  find  ten  algebraic  equations  to  solve  for  the  ten  unknov 
coefficients  in  the  zeroth  order  series  expansion,  two  sets  of  ten  algebraic 
equations  for  the  first  order  nearest  neighbor  modes,  etc. 

We  shall  illustrate  this  by  an  example. 
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III.  SPECIFIC  EXAMPLE  - CIRCULAR  CYLINDER 

A.  BOUNDARY  CONDITIONS  AND  EIGENFUNCTION  EXPANSIONS 

Let  us  consider  nonconcentr ic  circular  boundaries  such  that  (a,B)  = (r,0) 
and  (a1, S')  = (r’,9’)  with  r = a the  inner  and  r'=b  the  outer  boundaries. 


The  scale  factors  are  h^  = 1 and  h9  = r and  the  boundary  conditions  become 


(15) 


(16) 


(17) 


O = 


- ftoo1 


$ 4 


a_| 

dr 


A - A + A 1 
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r z ^11 
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2_  BAA  - 

T 1 B ev  rl  a© 


2. 

r 


d r 


& oor  [ 


A AT  ] 
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We  have  the  appropriate  eigenfunction  expansions 

(18>  I T.i.  Vr  ierA,..A«>A4*]fc"Ic«) 

<-)  IT,  [a.*. 


n = o 


<2°>  tf-sM  - ff'W  I [B,a.  T$».3,n.  3^]  {£]<«•»’> 

<21>  [ Y.  1 ,-.k  - Z "[  , Bwl  [-“I 


n=° 


which  reduce  the  boundary  conditions  to  algebraic  equations 

(22)  _ A x n.  - - f5  C O Ain.  + 2 M-  ^ A«4n  4 Ain.  4 or  Ash."  a.  Afett  ( 

(23)  O = 2ytL  £ - A^-  Asn.  “ 5Z  A-in  4 £ AfcU  J + 

(24)  _ ^ (VL  - - f5  CO1  S an.  + 2 M { “ ^ ®4n  * "b*  Bin.  ^ tZ  Bsn~  IZ  Bbn 

(25)  o = ? M [ ^ Bsn~'blBin.4”b  Bfen]  + fs  CO1-  Bsn 

(26)  5>2v\  r fV  G0L  4\r\  - "to  Bsn.1 
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We  now  consider  the  eigenfunction  expansion  of  the  Greens'  function: 

oo 

(27)  Q = - \]0  (4R)  = i Z COs[Tnfe-©o)]  Wm^p) 

4 0 m*  o 

2 2 2 

where  e = 1,  e =2  (m^O) , R = a +6  - 2aBcos(0-0  ) and  6 must  be  greater 

o m o 

than  a.  [These  8, a have  no  relationship  to  the  general  coordinates  used 
before  but  are  merely  convenient  labels  here  for  the  field  and  integration  points. 
0 or  0^  may  be  primed  or  unprimed  depending  on  the  particular  cases  and  a, 8 
may  be  either  a or  b . To  keep  a from  being  identically  equal  to  8 , we 
choose  field  points  just  "off"  the  integration  surface.  In  all  cases,  we  shall 
choose  them  to  be  "outside"  of  the  region  over  which  the  original  integration 
is  performed,  i.e.,  for  the  infinite  acoustic  fluid  we  choose  r'  at  b while 
for  the  elastic  solid  we  choose  r'  at  b+  and  r at  a-  for  field  points  at 
the  outer  and  inner  boundaries  respectively. 
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B.  BOUNDARY  INTEGRAL  EQUATION  - FLUID 

The  boundary  Integral  equation  for  the  infinite  acoustic  fluid  becomes 


£ yn  CoS  ( 0 - ) J 3"m  C k f ) Hm  (k  /o 

Yin-0 

(os [m  (e-0oj]  C v* ' ) b 

r=  b~ 
fo'=  t 


where  we  have  identified  a with  r’=b  as  smaller  in  magnitude  than  8 which 

is  identified  with  r ’=b.  This  tells  us  essentially  which  term  to  differen- 

o 

tiate  in  the  Greens’  function  series,  and  leads  to  a simple  relationship 

(29)  o=2  ? 3„(WlC(kUj(^J(«S'l 

Y)~0 

which  by  the  orthogonality  of  the  eigenfunctions  reduces,  for  Jn(kb)^0,  to 

(30)  Ba*  = k \-\a(M 


-1  u— 


C.  BOUNDARY  INTEGRAL  EQUATION  - SOLID 

W>  may  now  consider  the  boundary  integral  equations  for  the  elastic  shell. 
First,  we  examine  the  outer  boundary  equations  with  r'=b+,  i.e.,  we  place  the 
field  point  just  outside  of  the  outer  boundary  of  the  elastic  body.  Then,  for 


example,  we  have 


(31)  0 = 


2 _ 

s -* 


»>  7 b 6 ©o 

7 Jro'=b>v'=U- 


0 

s in 

The  first  integral  leads  directly  to  terms  involving  (cos?  n® ’ as  in  the 
acoustic  fluid  case  above.  The  second  integral  does  not.  We  cannot  mix  0 and 
0'  in  equation  (31)  as  long  as  they  have  different  origins.  To  avoid  this 
dilemma,  we  expand  about  6=0  such  that,  for  example, 

(32)  S &"rs=o]  S [H«<Fri]s=o 


The  S terms  may  be  removed  from  the  integration  leaving  an  integral  over 
C '*'(«-0),  i.e.,  some  Greens  function  with  0 and  0'  measured  from  the  same 


origin. 


Since  this  is  simply  a Taylor  Series  in  6 , we  have  (see  reference  [ ] for 


details) 


\ (V  Uo(fe^ 


(33)  G = j!  1 S S 


We  cannot  use  the  eigenfunction,  equation  (27),  expansion  of  G directly  since 
it  applies  only  when  0 and  0'  are  measured  from  the  same  origin.  Instead  we 


can  use,  for  example. 


(30  Js=, 


where 


r 


:35)  R*  = A + >ol-  2 A Ao  COS  ( 60-  © ) + 2 S [Ao  C05  00-  AtoS  © 3 + s 1 


[ S + Ao  cos  ©o  - A cos  © J 


Here  A and  are  distances  measured  along  8 and  0q  respectively  from 

origins  which  are  a distance  6 apart  on  the  x axis.  We  then  have  the  first 
order  kernel 

(37)  [ u*  0e1o]5=o  r - k [ V4. ( WR)  [Aoco5ec-Afoie]/R]5;o 


For  X - r , X = r f we  have 
o o 


08)  5 [H„  (**>]£’„  <W 

0 L 

TT  k An  (W)  [ 


We  may  set  rQ  equal  to  b and  r equal  to  a after  any  appropriate  differ- 
entiations are  carried  out. 

Similarly  for  A = r'  and  A = r 

o o 
z~n 

(39)  j [ [M0aR)];:o  d©, 

o 

. irlijnar.)[{s:]l'->*'l U„.,  M-[w] < 

Clearly  we  must  also  expand  all  of  the  other  dependent  variables  in  powers 
of  6 - or  equivalently  the  coefficients  of  their  eigenfunction  expansions: 


(a) 

A yn  ~ 

I 
1 '0 

4 . (A) 

A 

(b) 

00 

1 

i-o 

0 ^ f-O 

S 'Bjm 

-16- 


If  we  examine  the  boundary  conditions,  equations  (15-17)  and  the  acoustic 
fluid  boundary  integral  equation,  (20),  we  see  that  they  must  hold  for  each 
separate  order  (power  of  6)  since  5 does  not  appear  explicitly.  The  elastic 


boundary  integral  equation  (31)  however  becomes 


(41)  0 


(*"  f C r*  ? * 111  fsm? 

= ) S 2_  ^3rv  [(osj(neo) 

o j:<3  tv-o 

* 7 s‘Z  &’l  Ue°' 


Iso  rt;o 


r0'-b,lr'=f 


' £ s- £ A ;r  X s -s.-  ? ^ 

l — i — /-v  T - o , r = b 


1 = 0 v\  - ‘ 


which  may  be  expanded  in  powers  of  6 as 


0=  f fl 

o L n=o 

. ri[ 


r r!'(r[  n <■>  a&»  „ D ■"  fl  If  S'"  fine?)  7 
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2 r , <'l  i to  (of  A (01  9G-o"  A f0V  <"1/S"f,7(nfi.l 

♦ j l [A,.  £ - A<"  + A”  55-  -A«  G.  JfaMj 

° ■.xv0  . a d 9o 
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Solving  each  order  of  6 separately,  we  have 
(i)  Zeroth  order  equation 


(A3)  0= 


o r\-o 

(ii)  First  order  equation 

Oo 


<«>  HI  a«  e:“][s:l«"w 

o n=o 

- J " { | [-  *»«"  + &°]  l SJ  ("e;  > b de“  L.W.  r'.V 

, r r f [ a»"  - A«  6.“]{ -I'"*1  a4fcL-. ■-* 

o ^ n-° 

We  see  that  the  first  order  solution  is  coupled  to  the  zeroth  order  solution. 
Furthermore,  the  form  of  G is  such  that  an  j mode  zeroth  order  term 

will  produce  (j+1)  and  (j-1)  mode  first  order  terms.  A similar  equation  holds 

on  ip  . 

The  inner  boundary  integral  equations  are  handled  in  the  same  manner  with 
r’  equal  to  a , i.e.,  the  field  point  just  short  of  the  inner  boundary.  For 


example,  we  now  have 

( A5)  0 * 5 1 " 

O v. 


* Gc 


3$  (b) 

7 bdGo1 

arc' 

■7  r o' = b , r - a.' 

af(d)j 

f*  a d Go 

aro  J 

Vo=  0 > rr 

Now  the  first  integral  is  a mixed  type.  Using  the  same  expansions  as  before 
we  have 


A 
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(i)  Zeroth  order  equation 


o.  n f t-csf  • 

o^x\'o  r o - j ' - 

* ff  f [A,:’if:-A,rc.]fs} ^ 

f * — t-  r o ■"  o ^ v *■*  Q. 

O *»«»  ~ r* 


(ii)  First  order  equations 


(47)  0 = 


ff  z i-cf‘Cc  c 

O <f\  - 0 ro-'»irS® 

ffZK  ss'  *A\ 

J ' L r0^a>r=tL 


o y\ = 0 


with  a similar  equation  on  4’,  etc. 

We  now  choose  p(a)  to  represent  a single  mode,  i.e.,  A =0,  n i j 
and  A = P The  zeroth  order  boundary  integral  equations  and  boundary 

lj  .1 

conditions  will  also  only  involve  this  one  mode.  The  first  order  boundary 
integral  equations  will  be  satisfied  only  with  the  (j  + 1)  and  ( j — 1 ) modes; 
the  boundary  conditions  will  therefore  be  applied  only  for  these  modes  as  well. 
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D.  ZEROTH  ORDER  SOLUTION 

The  boundary  conditions,  equations  (15-17),  apply  for  each  order  separately 
and  may  be  used  directly  as  may  the  acoustic  fluid  boundary  integral  equation 
(30).  The  outer  boundary  shell  equation  on  $ , equation  (46),  integrates  to 

(48)  0=  I fee 

Since  H (kb)  cannot  be  zero  and  £sin{n0'  are  orthogonal,  we  have 
n D 1 cos) 

(»,  o . w { - b Kn  in  * a A,r  j/(k..)] ♦ t w (k» #) 

Similarly  for  ij,,  with  Hn(kRb)  not  equal  to  zero; 


(50)  0 


--  L {-1  BSr3/(h\»)+a  At*  Jn(kt) 


The  inner  boundary  shell  equations  give  similar  results  provided  J^dc^a)  and 


J (k„a)  are  not  zero, 
n R 


(50  o . W H C (k.M* o A Wo'  (H  * b B«  W.  lk.bl-0  AS  HnlM 
(52, 0.  - b c »; (wv).(x 
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In  summary,  our  equations  are 

(53)  (a)  -A«=  (-ps(0'  + 2*A«Va*)A»«^^C-A4v,/tt.^t,A5n/a*  -nAen/a] 

(b)  o = -vi/a-A  /a  4 /a.1  4 [&wV'a7o-1']  Asrt  4 A Afcn/a, 

z 

(c)  0=  4 [-&a>^?/uV/b\]BSn  •+?/u[-B-m/b4-nBSrv/bv-nB6^/b] 

(d)  0 = -nyu  B3n/b'  -*/M  BWb  4 C ^ Bsn  +/A  Bfcn/b 

(e)  0 = - B2n  v fF  60v  [ - Yl  Ben  /b  ] 

(f)  0 ~ " B art  4 & VW  (-kb)  8ta  / B\u  (bk) 

(g)  0=  b©  [-\j  ^Tn.  (if*  B 3n4  0-  Jyx  (ko  &)  Ain]  4 L Jr\  ( " Q. ZJ n (k©&)  A-tn 

(h)  Q = ()rr\>)  B sn  4 O-Tn  Asol  + ^ Xv\  (kR^  Bfin"  O-Tnlkaa)  Afcn 

(i)  0=  W[YYt  (kok)  Bsn  4 C^Yt  (k©0-)  Asnl  4 ^ Yn  (k*  k)  CL  Yn  U©a)  A«ri 

(j)  o-  Wf-VY/lWM  Bsn  4 a Yn  ( Ua-)  tUn]4k  Y,(feRk)  Bfc*  -a.  Yn(kRa)  Afcn 
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We  have  dropped  the  superscript  (0)  from  the  coefficients  in  these  equations 
and  trivially  combined  the  last  four  equations  to  give  real  rather  than 


complex  expressions,  although  Ain»  are  complex  anyway  due  to  the  complex 


term  in  the  sixth  equation.  If  A is  zero  for  n ^ j , all  coefficients  are 

In 


2ero  for  n ^ j.  This  then  leaves  only  one  mode,  n = j , to  be  determined. 
These  equations  may  readily  be  checked  by  the  well  known  solution  for  the 
concentric  case;  they  are  solved  here  to  indicate  the  procedure  to  be  used 
for  the  higher  order  terms  which  have  the  same  coefficients  but  a modified 
forcing  function  (inhomogeneous  terms). 


We  may  nondimensionalize  these  equations  using  a reference  "potential", 

2 2 

r*  1 t*  _ m.  _ /-\iO  \ i 


♦j  • 

related 

to  the 

reference 

We  use: 

(54) 

o(x>  - 

\?o  a J 

O^R  = B 

<2  = 

, - 

f ~ ) 1 

/(x+'zax')  - 2 


-i 


and 


(55) 


A,  » 

= o/t,  ■ 

B.  -*( 

B,-‘ 

Bzj  t>/ ( X+2/u) 

a5  = 

A»i/Si 

b3  = 

A+  - 

• *» 

• ro 

i-©i 

• •<+ 
< 

rw 

B^- 

■ b / $3 

As1 

A ^ j 

y 

A*,'- 

Bb= 

B>bJ  b / $ j 

leading  to 


J 
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(c)  o - e B,  4 ppx]  83  4 ^ [-  4 3 Bs  - j §<„"] 


(d>  o - - Bs  + Ba  4 [ p® " 3 B>s  4 ^ ® L 

(e)  o - - Bz.  -+  Pd  B4  " ^ pD  B>5 

/S/  ^ . 

(f)  0 = - Ba  -v 

(g)  0 = - p©  (po)  "Bi  + C^o3j  Wo)  As  4 J3  (jM  B4  - (°M  Aa 

(h)  0 = _p«  (^r)  B5  3^  W*0  B5  + (p«)  B>k  ~ WfO  At 

(i)  0 r ~ j3r>  Y^  ( p©)  Bs  * C^d  Yj  («/»)  As  4 ( pt>>)  B4  - Yj  (&t>)  A4 

(j)  o = -p*  Y/(pR)  Bs  4 *«Y/(°<*)As  + Yj(p*l  Bw-Y^  (^X«. 
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Although  this  may  appear  to  be  a formidable  system  - ten  equations  in 
ten  complex  unknowns  corresponding  to  a twenty  by  twenty  matrix  system  - 
the  equations  are  readily  reduced,  eliminating  variables  until  a simpler 
system  is  reached,  e.g.,  and  B,,  are  eliminated  directly  using  the  fifth 
and  sixth  equations: 

<■>  b*=  §5] 

(57) 

(b)  g,  = Pt  poL  £ $ * - 3 85] 

A3,  A5  can  be  written  in  terms  of  A4  and  A6  with  similar  equations  for 
the  B3,  B5 

A?  = + A->  + At 


(b) 

As  - 

Ks  + Csa  Aa-  + C5 1 Ab 

(58) 

(c) 

= 

DjA  E)4.  + T)3fc 

(d) 

Br  - 

D54  ,5a  + iPst  5 to 

using 

the 

first  four 

equations , 

V<j  = 

0/  C ( ~ X / C (O/o'-'t'-j 

with 

Ks  = 

1^»7c  - 1 

C|,  A 

= Cs  «,  = - 'fc'  cx't.V  C " 

1 

c 3fe 

= Cs^  - H’c1'5°,d1+i3 

V) /C  " 1 

and 

^34  = 

- p»  Cc-eSip»)/[(p»‘-tj‘) 

D 3 b : 

•*  a 


The  last  four  equations  (56  g- j ) are  rewritten  as 


~ ^ Q - - 


IT  o/o  io Alo  ^ 


Ao 


$5  ~ 


JlR 


As  = ■ gt  + At 

TTc^rSrvAr  o^r. 

Ap((X,M=  ( a)  Y^  (b)  - 3~j  ( b)  Y j T<X) 

( a , b}  = J^(cO  ( b)  — (b^>  Y}  Ca) 

Jl  c<x,  b)  = Xj  7(  a)  (b)  - Tj  (bl  Y3  Co.) 

A.(a>M=  Tj  (a-)  Yj  (b)  - XjYb)  Y^  ('a-') 


where 


and  pp  = pCa^.Bp)  etc.  We  then  have  a four  by  four  system.  Using 
Kp  = 2/  (wotpBpSp)  and  = 2/ (ttcir0rSr)  we  have. 


0 D34I  ’^3C  A*< 

0 -1^54  [^R  Ab 

[ A*/*,  A.  - Y "Cw  ■ Ko  0 B‘ 

-Css  [ Csfc]  0 -Xr 


with  a solution  given  explicitly  in  appendix  A. 
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By  backwards  substitution,  we  can  evaluate  all  of  the  unknown  coefficients. 
These  must,  of  course,  agree  with  the  well-known  solution  for  the  concentric 
case  and  are  carried  out  in  detail  to  serve  as  a procedural  check  for  the 
higher  order  solutions. 
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E.  FIRST  ORDER  SOLUTION 

Again,  we  use  the  boundary  conditions  and  acoustic  fluid  boundary  integral 
equation  directly  for  the  first  order  coefficients.  The  outer  boundary  shell 
integral  equation  on  , equation  (44),  integrates  to 

oo  t- 

o = Z ~ B 3*  ^ J**>  ' 2TY  3n.  k 3V  (kiA)  j \-\r\  (ioU 

(62)  „.'o  J 

+ | a ko  (k11)  " 0-  ?-TY  Tn  ( VW  (feok) 


Z 0-  TV  W 3yi  0?o  ft-)  - Aa*  a IT  kt>  (1*0®)  j 

n^o  I L 


and  a similar  equation  on  i|j  with  A^,  A^,  B^,  B^  and  kR  replacing 

A3n’  \n’  B3n’  ®4n  and  S)  * 

The  inner  boundary  shell  integral  equation  on  <J>  becomes 

(63)  0*  Z lf“  ^ \>  ko  'Tf  Wy\  0?  ok)  + B«>n  1 V ko  VI « ( kc 

r>=o  l 

• { X,.,  (!..«)-  { T.-, »-«)]  ] + 


* Z W A*  Wr\  + ^ 2VY 

+ C'ft  ^ U Wn(U-A^'  ft  ZTt  VA«Uo<o] 

• [ 3r>(UaV  V'Jl]^e'} 


with  again  a similar  equation  on  <p  . 

Clearly  if  there  is  only  one  zeroth  order  mode,  n = j , there  are  only 
two  first  order  modes,  n = j+1  and  n = j-1  . Consider  the  j+1  mode  first. 
The  governing  equations  are 


/ r>) 

63  ^+1  k’o^1  3"^.+. , C V?t>  ^3  3 ■+•  -j*  > ' b ■ 3”^  + | ( ko  1 

. (<>  . __  / (O 

/''3>3+l  Rp  CL ' J j+i  (1*,>Q.)  — 4j  + i ' (X  3j+i  (ited) 

2 I ^ !?» Q-  3^ 

k [ B3'°’  Ut  j/(l..b3-  Bi^’  Lb  J,(hl>)] 


CO  t t * ">  , , . , . . 

“Bb^+i  \?b  t V~\  ( to  b ) * 3 4 + . ( to  t» ) 

+ A^+»  to  a V\  jti  (to  <3  - Aa^<  a ^»a) 

i [ s5f  W.‘W-  U/(WM-  B+'"'  UL  Hj  (U)] 
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and  two  similar  equations  on  0 . The  j-1  mode  is  essentially  the  same  with 
a sign  change  for  the  inhomogeneous  (zeroth  order)  terms 

(66)  - 8^-,  WV>Tj_,  w 

+ Wo  a J (Wo©-)  - Aa  01 ' 

- X [ Avj  ^ Wo*  a A43  ’ Wo  a (kDa)] 


— £>3^-1  • Wo W ■ \-\ j-t  (WP\) ) + B-, 3_,  W M^-i(WoW) 

♦ AiV-1  ^V»(koQ-)  - A44-1  a H -3 -•  ( Wo <^) 

- T [ B330'  WD*W  ^^(WoW)-B43^0)  tot-  Uj  ( WoW)"] 

and  two  similar  equations  on  t|i  . 

We  note  that  the  homogeneous  part  of  these  equations  plus  the  first  order 
boundary  conditions  and  the  acoustic  boundary  integral  equation  are  identical 
to  the  zeroth  order  equations  with  j+1  and  j-1  for  the  first  order  coeffi- 
cients replacing  j for  the  zeroth  order  coefficients.  The  inhomogeneous  terms 
are  modified  and  now  appear  only  in  the  four  shell  boundary  integral  equations. 
Then  we  can  reuse  the  solution  scheme  for  the  zeroth  order  j mode  coefficients 
to  solve  for  the  first  order  (j+1)  and  (j-1)  mode  coefficients  with  a minor 
change.  We  again  may  write  equations  (65)  and  (67)  with  real  coefficients  by 
combining  with  the  previous  equations  - equivalently,  we  replace  Hn  by  Yn 


in  these  equations. 
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We  may  again  nondimensionalize  using  an  additional  length  scale  factor  of 
a for  these  first  order  terms  which  must,  in  final  form,  be  multiplied  by  the 
length  6 in  order  to  be  compared  to  the  zeroth  order  terms;  e.g.. 


<Ay  ( 1 ) 

As54, 

A 

- 

*/*i 

^ ( ' ) 
A-4  ^41 

A 

- A-4-J4I 

( > ^ 

- 

a/i-i 

B4 

( » > 

= 84.^^ 

cr 

N. 

et- c ■ 

scale  factor 

is  used  for 

both  A(1) 

of  equations  as  in  the  zeroth  order  case.  We  then  have 

, ~ C>»  r ~ r (i)  ~K  Ol 

(a)  0 - [ t ( ^ + ')  * ] Al^>  «■  t + ( As  J 

(b)  0 s-r  ( 3+0  Avj+.+rl^')  A^-m  + [ - (j-h )r]  As  *-m  ♦ ^ At  *♦» 

(c)  0=  ~ ^ + ( Bs  j-n  “ ( 6b.j4i  ) 

(68) 

(d)  o=  - ♦[ 

Ay  m\  /.i  ^ 

(e)  o = - -*■  [ 8<v$*.  -(^+1)  Bs^+i"l 

A^  (I)  ^ ( l>  / 

(f)  o - - 8M+.  + 

(g)  Li  = ^ [ p®  (fM  8-^  - 

~ <.)  _ / 7 «»  T / 7 «•»  -r  x 7 «•>  T , 

* “*  D 3^+t  J^c  (<Vo^  * U4j4»  (°^1 


1 


•30 


(h) 


V.,  = ^ [ 8*7  J3Vl  - Bb'j' 

” ~ PR  + by*  °^R  Y^+i  St, j+t _3^+l  (pR) - At,^+i"3^i  lot*) 

(i)  L*  = T \ *>}$  |3t>Y^(^o)~  B+j  Y^  (fii^l 

r_  P°  Y-J-**  ■*■  ^iy'  Y^,  (°/t>)  * 84  Y}+,  (jio)'  Aajv\  \^4v  (o/d) 

(J)  \_  a - °r  T Bsj  PR  Yj  ( p*)  - Bt^  Y^  ( pp')"] 

= - Bs^»  p«  Y^,  (jin)  + ASi+,  c/*  Y^v>  (°(f0  + Bt  Y^+i  (jjrV  Afcj+iY^ti  (°^/ 


and  a similar  set  of  equations  with  j-1  replacing  j+1  and  -L^,  -L^.  -L3’ 

-L4  replacing  L , L2>  L3>  L^. 

We  see  at  this  point  that  we  have  two  separate  ten  by  ten  systems  of 
complex  algebraic  equations  - each  very  similar  to  those  of  the  zeroth  order 
case.  We  again  may  reduce  the  unknowns  to  yield  a feasible  arithmetic  problem  - 


consider 

the  j+1  system 

as  an  example: 

(a) 

‘X  ('•> 

= p 

: (s:;:.  - is*" 

§ 

-7 

5-J+'  J 

(b) 

B » £*>  ~ 

[ 

- 

"I 

Bsj-H  J 

(c) 

r i" 

A3^  + . - 

c5* 

A4^+. 

-4- 

C^fc 

A 

AtJ4, 

(69) 

(d) 

A 5$+'  - 

C 54 

C < ' 

A<i>' 

+ 

C 5 fc 

A«  -J4( 

(e) 

Q - 

°3^'  “ 

Y>>a 

e^v. 

+ 

Dik 

Q 

Dfe 

(f) 

<N/  . . 

B,V-  * 

_ X <«  i 

D54 

-v 

T)  s c. 

~ t<  ) 

Bfc^» 
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where  the  C's  and  D's  are  the  same  as  before  with  j+1  replacing  n . 
Clearly  the  and  K,.  are  no  longer  present  since  the  inhomogeneity  in  the 

system  has  moved  to  the  last  four  equations. 

Those  last  equations  are  now 


(a)  VU  Via)  ~ D3  fc  6b^.  = 

~ ..x  _ ^ o'  I As  \ r?  O'  Xjai  i (<*«!  _ p 

hi  8< * (iT.,-  M 6t)..  = —J-^TZ- Kl 


(b)  Vis  Abj 


(70) 


\ v _r  ir"  r a“'  w r(°  - p 

(C>  lo/oM  ^a]  A^4,  -UfaAb^'  " V(o  044*1  " 


/ jn  _ P A11’  * r 5*  - r I a “ 1 -\/  Q<0  - La  L»  Yyi  (p«)_  p 

(d)  Ls*  A44A'  UtJ  Abj*'  ^ 


remembering  that  the  C,  D,  p,  q,  r and  s are  all  modified  to  have  j+1 
replace  n . 

The  equations  for  the  (j-1)  mode  are  identical  to  the  above  with  the 
corresponding  changes  of  j-1  for  j+1  everywhere  and  with  the  sign  of  the 
R terms  reversed.  The  solution  procedure  is  the  same  as  for  the  zeroth 
order,  and  is  also  given  in  appendix  A. 

To  convert  these  first  order  coefficients  to  a form  suitable  for  compari- 
son to  the  zeroth  order  solution,  they  must  be  multiplied  by  (6/a),  which  is 
a small  number  in  order  to  have  the  first  order  be  a sufficient  correction  to 
the  zeroth  order  solution.  Higher  orders  may  be  calculated  in  the  same 
manner  — only  a new  set  of  Lj  need  be  calculated  since  the  same  solution 
procedure  will  hold  for  higher  orders  as  did  for  the  first  order,  i.e.,  the 


same  basic  program. 
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CALCULATIONS  AND  CONCLUSIONS 

The  primary  application  of  this  solution  is  to  a thin  steel  shell 

submerged  in  water.  The  parameters  used  were  a = 300  cm,  b = 303  cm, 

c = 146000  cm/sec,  c = 496000  cm/sec  and  c = 270000  cm/sec  with 

I)  K 

3 3 

P = 7.7  gm/cm  and  p = 1.0  gm/cm  . Results  for  the  exterior  pressure  B , 
Sr  1 

are  shown  in  Figures  2 and  3 for  the  case  j = 0 and  j = 2 respectively  for 
a unit  interior  forcing  pressure.  Results  for  the  first  order  solutions 
are  given  on  the  same  figure  as  the  zeroth  order  solution,  but  to  a 
scale  50  times  as  great.  Since  the  first  order  terms  must  be  multiplied 
by  6/a  to  determine  their  contribution  to  the  exterior  pressure  field, 

_3 

the  actual  exterior  pressure  will  depend  on  6.  For  6 = 0.3  cm,  6/a  = 10 
and  the  first  order  contribution  is  generally  negligible.  There  are 
regions,  however,  where  the  zeroth  order  solution  is  minimum,  e.g. 
near  300  rad/sec  for  j = 2.  In  this  region,  the  first  order  contributions 
may  be  comparable  to  the  zeroth  order,  but  this  is  for  a very  limited 
frequency  band.  Figure  4 shows  results  for  a thick  shell  with  a = 200  cm 
and  the  remaining  parameters  unchanged.  In  all  cases,  w runs  from  200 
to  9000  rad/sec.  Here  the  first  order  and  zeroth  order  solutions  are 
plotted  on  the  same  scale  and  are  very  similar.  This  implies  the  first 
order  contribution,  when  multiplied  by  6/a  will  be  negligible  for  any 
example  in  which  6/a  is  small  enough  to  allow  a first  order  theory  to 
be  used. 

The  general  conclusions  to  be  drawn  from  this  work  are  first  that 
such  an  approach  is  applicable  to  radiating  submerged  "nonconcentric" 


elastic  shells  and  second  that  nearest  neighbor  modes  are  excited  in  the 
fluid  in  the  higher  order  solutions  for  a single  mode  forcing  function. 
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Both  of  these  conclusions  have  already  been  drawn  in  the  acoustic  shell 
case,  reference  [1];  we  have  shown  here  that  they  are  applicable  to  the 
physically  more  meaningful  elastic  shell  case  as  well,  with  some  addi- 
tional effort.  A specific  result  based  on  these  two  examples  is  that 
the  influence  of  nonconcentricity  is  essentially  negligible  in  most  cases 
where  6/a  is  small.  Since  higher  order  solutions  may  be  computed  with 
the  same  basic  program,  extensions  to  larger  values  of  (6/a)  are  clearly 
possible  for  thicker  shells.  Such  analytical  solutions  may  be  useful 
both  in  their  own  right  and  even  more  importantly  as  checks  on  general 
numerical  schemes  which  might  not  be  otherwise  properly  constructed  to 
include  (or  recognize)  such  effects  as  nearest  neighbor  modes. 
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APPENDIX  A 
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A4j  = N4/D 
^(o) 

A6j  = N6/D 
(First  order) 


D (j+1)  = D with  j+1  replacing  j in  all  terms 


N 4 (j+1): 

R^,  R2 . R4  defined  in  equation  70 

N6  (j+D: 

and  a similar  set  for  j-1  contribution. 
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A3j 

- k3  + 
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A5j 
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CS4  A^j  + C56  A6j 
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% (0) 
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APPENDIX  B 


A simple  test  calculation  is  the  "transparent"  shell,  i.e.  pp-  Pg, 
= o and  cQ  = c.  Since  the  pressure  is  applied  at  the  inner  surface. 


r = a,  we  have 


in  V\n"  (kr)  f sin) 

Vr>*]  = — 


On  an  "offset"  circle,  we  have  r'  = b given  by  r = r ' + 6 cos  6 +. . . 


r = b + 6cos  6 . 


e)=  f h ^ cos&tvCO^Vb  ms  J1 

V-C’ftO  L 


where 


\AvT0^  . 


% (*,&>  = i>n 


V (VI  -- 


If  we  return  to  the  integral  equation  approach,  a direct  substitution 


of  these  parameters  into  the  zeroth  order  equations  leads  to 


O "»  » Hsifil 

V\y>  (°0 


I 
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The  first  order  equations  lead  to 

(l)  V 1HI  P 

ft'nu  = '1  vC’  (V*)  " 


Q _ _* 

D ' »-  ‘ z 


* \W  t'"' 

1 H.1-',’  (M3  Jn 

2 W"(M> 


We  must  bear  in  mind,  however,  that  these  are  values  for  the  pressure 
at  the  outer  surface  measured  in  terms  of  0'  and  not  in  terms  of  0.  We 


need  to  convert 


(sln7  (n 
I COS  \ v 


0.),  We  need  to  convert 


{£:?(»  e)to[^ e') which 


modifies  our  first  order  results.  We  have  the  definition 


01-  0+  sm"*  £ ^ s\y\oJ 


and  to  a first  order 


{ 22]  = {23]  <"•'»  - * Ph>6]i 


b:;’ 


" L t0SJ  s w 


VU,  M 


* * vu  w ♦ = 


which  agrees  with  the  previous  solution. 
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